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ABSTRACT 

HV  used  the  recently  developed  S-tmusform,  an 
extenfiryn  to  the  ideas  of  the  Gabor  inmsfonn 
and  wavelet  transform,  to  analyze  the  Tropical 
Ocean  Global  Atmosphere  (TOGA)  sea  level 
data  set  und  to  obtain  the  localized  spectra, 
varying  with  time.  Our  S-spectra  show  some 
features  tha  t  can  ant  be  obtained  from  the  Fou¬ 
rier  transform ,  such  as  phase  lock,  temporal 
reliability  of  spectra,  and  out-of-phase  behav¬ 
ior  of  the  semiannual ,  annual,  and  biannual 
signals  between  the  western  Pacific  (Maura) 
station  and  the  eastern  Pacific  (1m  Libertad) 
station-  The  annual  signal  is  usually  pliasc- 
lockcd  with  the  semiannual  signal,  but  not 
with  tin-  biannual  signal.  The  annual  signal  is 
guile  stationary  in  the  western  Pacific  station 
and  nan-stationary  in  the  eastern  Pacific  sta¬ 
tion.  Du  ring  1980-84.  the  quasi-biannual 
( QB)  signal  was  very  strong  at  the  western 
Pacific  station  and  quite  weak  at  the  eastern 
Pacific  station  However,  during  1974-70  and 
1986-90.  the  QB  signal  was  weak  at.  the  west¬ 
ern  Pacific  station  and  strong  at  the  eastern 
Pacific  station  This  may  imply  different  phys¬ 
ical  processes  involved  in  the  western  and  (ast¬ 
ern  Pacific  during  the  El  Nino  and  Southern 
Oscillation  lENSO)  jMfriods. 


INTRODUCTION 


After  Pedlosky  ( 1970 )  proposed  the  dual  time 
scale  concept  for  nonlinear  waves,  a  one¬ 
dimensional  perturbation  along  x-axis  is  gener¬ 
ally  depicted  as 


tp  =  A(io,  t)  C0S( k.v  —  2n-jt.)  (1) 

where  k  is  the  wave  number,  u>  the  frequency, 
and  (x.  t)  the  dual  time  scales.  The  temporal 
variable  t  represents  the  short  time  scale  (fast 
variation),  and  t  denotes  the  long  time  scale 
(slow  variation).  .4  is  the  amplitude  which  var¬ 
ies  on  a  slower  tune  settle  r. 

Many  papers  have  been  published  since 
then  to  determine  ,4(<a,  7)  from  different  dynam¬ 
ical  models,  especially  from  finite-amplitude 
wave  models  (e.g..  Pierrehumbert  1984;  Ped- 
loskv  1992).  However,  little  work  has  been  done 
to  obtain  A( 01,  7)  from  real  ocean  data.  The 
FVturier  transform  is  more  commonly  used  to 
obtain  the  spectrum  from  ocean  data.  The  mutual 
relation  of  the  Fourier  spectrum  H(tn)  and  its 
time  series  h(l)  is  given  by 

Ufa)  =  f  It  (if  ^dt  (2) 


and 


h(l)  =  ^  J  H(tn)eciwmid<u  (3) 

Obviously,  the  Fourier  spectrum  Ilf  01)  does  not 
depend  on  time.  Statistically  the  amplitude 
A( to.  7)  in  (1)  is  considered  as  the  absolute  value 
of  a  localizing  spectrum  near  time  t.  To  analyze 
the  time  variation  of  the  sea  surface  height  spec¬ 
tra.  we  need  to  use  some  other  transform,  such 
as  a  wavelet  transform,  the  Si  tort  Time  Fourier 
transform  (ie.,  the  Gabor  transform),  or  the 
S-transform. 

Tlie  Gabor  transform  is  given  by 
(Gabor  1946) 


I  You,  7)  =  J  hft)g(t  -  7.  <r)e  '-wrff  (4) 

where  g(t,  aj  is  called  the  Gaussian  Window, 
defined  bv 


<]( t,  a) 


1 


V  2xr*T 


exp 


-(-u) 

\V2PJ 


(5) 


where  x  is  the  translation  parameter  with  the 
same  dimension  as  I.  and  it  is  the  window  width. 
For  any  x.  the  Gabor  transform  provides  tite  spec¬ 
trum  of  a  windowed  segment  of  hft),  centered 
at  time  x.  The  Gabor  transform  has  two  kinds 
of  limitations.  First,  if  the  dataset  has  a  transient 
component  with  a  scale  smaller  than  it.  it,  is 
difficult  to  locate  it  with  precision  better  titan 
it.  Second,  if  the  data  set  has  important  features 
1  »f  differing  sizes  then  an  optimal  g<  l.tr)  for  analy¬ 
sis  cannot  be  found.  Therefore,  the  Gabor  tixuts- 
form  is  more  suitable  for  analyzing  data  where 
all  features  appear  approximately  at  the  same 
scale. 

The  wavelet  transform  W(~.  a) 
addresses  the  resolution  problem  by  introduc¬ 
ing  a  dilation  (or  scale!  parameter  a  and  is 
defined  by 


W( 7,  a)  -  J  h(t)w(t  —  x,  erjdl  (6) 

where  the  functions  w(t ,  a)  are  called  wavelets. 
The  spectrum  of  h(l)  Ls  obtained  througli  corre¬ 
lation  or  convolution  with  w(t,  a).  The  dilution 
scale  it  determines  the  width  of  the  wavelet 
w(l,  a)  and  thus  controls  the  resolution.  In  the 
wavelet  transform,  the  frequency  is  not  explic¬ 
itly  presented.  There  is  no  direct  relationship 
between  Fourier  spectrum  and  wavelet  spec¬ 
trum. 

Different  from  the  Gabor  transform 
and  wavelet  transform,  the  S-transform  lias 
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variable  windows  and  a  close  connection  to  the 
Founer  transform.  In  this  paper,  wc  use  the 
S-transform  to  obtain  a  localized  spectrum  of 
the  sea-level  front  the  Tropical  Ocean  Global 
Atmosphere  (TOGA)  Data  Set. 

THE  S-TRANSFORM 

he  S-transform,  recently  proposer!  by  Stock- 
well  et  al.  (1994).  is  a  generalization  of  the 
Gabor  transform,  and  an  extension  of  the  Wave¬ 
let  transform.  The  mutual  relation  of  the  S-spec- 
trum  S(o>,  t)  and  its  time  series  h(l) is  given  by 

SYoi,  t)  =  J  Il(u>  +  ajr  ^  ?'-™da  (7) 

and 


A, ft)  =  sin(lOf)  +  sin(20f).  0  £  t  ^  200 
and 

h/t)  =  2sin(20t),  0  £  f  <  100 
fcj(100)  =  sin(2000)  -  sin(1000) 
h/l.)  =  2sin(10/),  100  <  t  s  200 

The  first  signal  consists  of  superposition  of  two 
frequencies  (Figure  la),  and  the  second  signal 
consists  of  the  same  two  frequencies,  each 
appearing  separately  over  half  of  live  signal  dura¬ 
tion  (  Figure  lb).  The  magnitudes  of  the  two 
Fourier  sjtectra  are  identical  (Figure  lc.  d).  indi¬ 
cating  the  incapability  of  distinguishing  the  two 
signals.  On  the  other  hand,  the  magnitude  of  the 
two  S-speetra  are  very  different  (Figure  le.  f). 

TOGA  SEA  LEVEL  DATA 


h(i)  =  J  J  S(  o>,  rj^OrdM  (8) 

where  H(u>  -  a)  is  the  frequency-shifted  Fourier 
spectrum  of  hit).  Comparing  Equation  (S)  with 
(3)  one  can  lind  that  the  Fourier  spectrum  H(  u>) 
is  a  time  average  of  the  S-spectrum  ,S'( u>,  -): 

H(u>)  -  J  .*?(<■>.  t)(/t  (9) 


Equation  (8)  can  also  be  viewed  as  the  decompo¬ 
sition  of  a  time  series  h(t)  into  sinusoidal  oscilla¬ 
tions  on  time  I  with  the  temporal  varying  ampli¬ 
tudes  .SV  io,  "). 

Let  h(k&T ),  k  =  0.  I . A’  -  1  denote 

time  series,  corresponding  to  h(t),  with  a  time 
sampling  interval  of  AT.  The  discrete  form  of 
Equation  (7)  is  given  by 


The  monthly  sea  level  data  provided  by  the 
TOGA  Sea  Level  Cenrer  at  the  University 
of  Hawaii  is  list'd  to  demonstrate  the  usefulness 
of  the  S-transform.  For  simplicity,  only  two  sta¬ 
tions  for  the  computation  are  selected:  Nauru 
1974-90  (0'32’S.  166u.>U£r)  and  la  libertad 
1971-89  (2C12  S,  80*55' H)  representing  equato¬ 
rial  western  and  eastern  Pacific  (Figure  2). 

There  are  some  missing  observations 
in  the  two  stations  (Figure  3).  Complete  time 
series  (no-missing  data)  for  computing  the  spec¬ 
trum  is  needed.  Therefore,  the  linear  interpola¬ 
tion  to  obtain  values  for  the  missing  data  is  used. 
From  the  modified  time  series  of  sea-level  Fig¬ 
ure  -1),  annual  and  interannual  variabilities  art- 
evident  at  both  stations. 


FOURIER  SPECTRA  OF  SEA 
LEVEL  DATA 


The  sea  level  data  have  been  detrended  and 
nondimensionalized  by 


for  n  4-  0  and  by 

•V-  I 

«> 

for  n  =  0.  Here  the  discrete  Fourier  spectrum  is 
computed  by  (Brigham  1974): 

V-l 

jaUrl  =  T-SfdkAT)^  (12) 

[A- A  I  A  m  r* 

LIMITATION  OF  THE  FOURIER 
TRANSFORM 

he  spectrum  computed  from  the  Fourier 
transform  only  shows  overall  beltavior  of 
a  data  set  Different  data  sets  might  itave  the 
same  Fourier  transform.  For  example,  consider 
the  following  two  signals: 


h  -  li 


before  computing  the  spectra.  Here  h  is  the 
mean  value  of  the  sea-level  at  the  station.  Due 
to  the  short  time  series,  the  Fourier  spectra 
( Figure  5),  are  aseless  for  periods  beyond 

sixty  months.  All  the  interesting  peaks  are 
crammed  into  the  domain  with  periods  shorter 
titan  thirty  months:  (1)  several  signals  (semi¬ 
annual.  annual,  and  bi-annual )  appearing  in  both 
Nauru  and  La  Libertad  spectra:  (2)  a  strong 
annual  signal  in  the  western  Pacific;  and  (3)  a 
strong  bi-annual  signal  in  the  eastern  Pacific. 
Titus,  the  Fourier  spectra  provides  overall  infor¬ 
mation  (in  the  sense  of  time  averaged  spec¬ 
trum)  about  the  sea  level  variability.  Do  these 
signals  (semi-annual,  annual,  and  bi-annual) 
Itave  constant  amplitudes  during  the  whole  sam¬ 
pling  period?  To  answer  this  question,  we 
should  use  the  S-transform. 
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magnitude  of  tfie  Fourier  si 
e  h2. 


S-SPECTRA  OF  SEA  LEVEL  DATA 


Using  Equations  (10)  and  (11)  we  compute 
the  S-spectra  from  the  nondimensional  sea 
level  data  for  the  two  stations  i.  Figure  2).  S  is  a 
complex  function  of  ut(  =  wArAT)  and  time 
-(  =  j\T).  Here  A  7"  =  1  month.  The  period  P  is 
defined  by 


(13) 


If  we  use  the  unit  of  year  for  period.  Equation 
(13)  becomes 


(14) 


Thus  the  complex  function  S  can  be  transformed 
into  a  function  of  P  and  t.  The  amplitude  and 
the  phase  of  S(P.  t)  are  computed  by 


Figure  3.  Original  monthly  sea  level  (mm)  at  (a)  Nauru  (western  Pacific),  and  (b)  La 
libertad  (eastern  Pacific). 


A(P.  T)  =  1 SCP,  1)\, 


<b(P.  r)  =  arctan 


hn[S(P.  r)l 
H'iSfP.  7)}. 


(15) 


The  function  A.  called  the  S-spectrunt.  is  plotted 
against  the  (log.  P)  and  time  t  for  Nauru  and 
La  Libertad,  as  shown  in  Figure  6.  If  we  average 
A(P.  7)  along  die  --axis,  we  will  get  the  Fourier 
spectra  Figure  6  shows  several  interesting  fea¬ 
tures  at  both  stations:  (1)  the  annual  signal  was 
phase  locked  with  the  semiannual  signal  in  both 
the  western  and  eastern  Pacific  since  the  maxi¬ 
mum  values  of  A(  1,  t)  usually  correspond  to  the 
maximum  values  of  .4(1/2.-):  and  (2)  the  phase 
lock  of  the  annual  and  bi-annual  signals  is  not 
very  evident  except  during  1974-76  in  the  wesl- 
em  Pacific  (Figure  6a)  where  high  values  of 
A(P,  7)  extend  from  P  =  1/2  year  to  P  =  1  year. 
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Figure  6.  S-specira  normalized  by  the  maximum  value  at  (a)  Nauru  (maximum  value  5.59).  and  (b)  La  Libertad 
(maximum  value  =  2.35). 
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Figure  7.  Voices  of  the  semiannual-  oscillation  for  (a)  Nauru,  and  (fa)  La  Libertad. 


SEMIANNUAL  AND  ANNUAL 
SIGNALS 


Since  A(P,  r)  and  4>fP,  ~)  are  the  amplitude 
and  phase  of  a  complex  variable,  St  P.  t).  the 
sinusoidal  function 


V(P.  7)  =  A(P.  t )* 
COS  I  ^  +  I\)(P.  7) 


(16) 


provides  quite  useful  information  about  the  time 
evolution  of  the  S-spectrum.  The  function. 

V(P,  7),  evaluated  at  a  particular  period  P.  is 
called  the  “voice."  The  voice  for  P  =  1/2  year 
(semi-annual  signals)  shows  some  interesting 
features  at  the  western  and  eastern  Pacific  (Fig¬ 
ure  7):  (1)  the  weakest  semiannual  signal 
is  found  in  1977  in  the  western  Pacific  and  in 
1981  in  the  eastern  Pacific;  (2)  the  strongest 
semiannual  signal  is  found  in  1984  in  both  the 
western  and  eastern  Pacific;  and  (3)  during 
1974-80  the  strength  of  the  semiannual  signal 
varies  rapidly  in  the  western  Pacific,  but  is 
quasi-slationary  in  the  eastern  Pacific. 

The  voice  for  P  -  1  year  (annual  signals) 
indicates  that  (Figure  8):  (1)  the  annual  signal 
is  quite  stationary  in  the  western  Pacific  and 
non-stationary  in  the  eastern  Pacific;  and  (2) 
during  1971-74  and  1981-86,  the  annual  signal  is 
very  weak  in  the  eastern  Pacific. 


(a) 


(b) 


QUASI-BLANNUAL  (QB)  SIGNALS 

When  the  voice  is  evaluated  at  particular 
frequencies,  uh  =  S/A'AT.  the  function 
l‘( ox.  7)  represents  QB  signals:  25.5  month 
period  at  Nauru  (Figure  9a)  and  28.5  month 
period  at  La  Libertad  (Figure  9b).  The  QB  signal 
has  an  out-of-phase  feature  at  the  western  and 
eastern  Pacific;  large  amplitudes  appear  during 
1980-84  at  Nauru  and  during  1972-73  and  1986- 
89  at  La  Libertad.  and  small  amplitudes  emerge 
during  1974-77  and  1987-90  at  Nauru  and  during 
1977-82  at  La  Libertad.  This  may  imply  different 
physical  processes  involved  in  the  western  and 
eastern  Pacific  during  the  El  Nino  and  Southern 
Oscillation  (ENSO)  periods. 


CONCLUSIONS 

(1)  Our  S-spectra  show  some  features 
thai  cannot  be  obtained  from  the  Fourier  trans¬ 
form,  such  as  phase  lock,  temporal  variability 
of  spectra,  and  out-of-phase  behavior  of  the 
semiannual,  annual,  and  biannual  signals 
between  the  western  and  eastern  Pacific  stations. 

(2  )  The  annual  signal  is  usually  plia.se- 
locked  with  the  semiannual  signal,  but  not  with 
the  biannual  signal.  The  annual  signal  is  quite 
stationary  in  the  western  Pacific  station  and 


MTS  Journal  •  VoL  29,  No.  4  •  35 


Voice  (n=8) 


non-stationary  in  the  eastern  Pacific  station. 

(3)  We  found  out-of-phase  behavior  of 
the  QB  signal  in  the  western  and  eastern  Pacific 
stations.  During  1980-84,  the  QB  signal  is  very 
strong  at  the  western  Pacific  station  and  quite 
weak  at  the  eastern  Pacific  station.  However, 
during  1974-76  and  1986-90,  the  QB  signal  is 
weak  at  the  western  Pacific  station  and  strong 
at  the  eastern  Pacific  station.  This  may  imply 
different  physical  processes  involved  in  the 
western  and  eastern  Pacific  during  the  ENSO 
periods. 

(4)  The  statistical  significance  and  con¬ 
fidence  interval  of  the  S-spectrum  will  he  dis¬ 
cussed  in  another  paper. 
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